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Abstract 

The optimization of the usual entropy Si\p] = — J du p{u) In p{u) un- 
der appropriate constraints is closely related to the Gaussian form of the 
exact time-dependent solution of the Fokker-Planck equation describing an 
important class of normal diffusions. We show here that the optimization 
of the generalized entropic form Sq[p] = {1 — / du [p{u)]''}/{q — 1) (with 
q = 1 + — 1/ G 7?.) is closely related to the calculation of the exact time- 
dependent solutions of a generalized, nonlinear, Fokker Planck equation, 
namely ^p^ = ~-§^[F'ix)p^] + ^^^P^ ^ associated with anomalous diffusion 
in the presence of the external force F{x) = ki — k2X. Consequently, paradig- 
matic types of normal {q = 1) and anomalous {q ^ 1) diffusions occurring in 
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a great variety of physical situations become unified in a single picture. 
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Anomalous diffusion is intensively studied nowadays, both theoretically and experimen- 
tally. It is observed, for instance, in CTAB micelles dissolved in salted water the analysis 
of heartbeat histograms in a healthy individual 0, chaotic transport in laminar fluid flow of 
a water-glycerol mixture in a rapidly rotating annulus [Q , subrecoil laser cooling |Q , particle 
chaotic dynamics along the stochastic web associated with a d=3 Hamiltonian flow with 
hexagonal symmetry in a plane f^, conservative motion in a d=2 periodic potential [§, 
transport of fluid in porous media (see and references therein), surface growth 0, and 
many other interesting physical systems. Its thermostatistical foundation (as it is known 
for normal diffusion) is naturally highly desirable and has, since long, been looked for (see, 
for instance, and references therein). This goal was recently achieved by Alemany and 
Zanette (see also [1^) for Levy-like anomalous diffusion, in the context of a generalized, 



not necessarily extensive (additive), thermostatistics that has been recently proposed |1TT| , [l2 



This thermostatistics (already applied to a considerable variety of physical systems [jl3[ and 
optimization techniques |0) is based upon the entropic form 

which reduces, in the g — > 1 limit, to the standard Boltzmann-Gibbs entropy 

Si[p] = ~ J du p{u) Inpiu) (2) 

We show here that an approach similar to that of |^ makes possible a quite satisfactory 
discussion of a sensibly different anomalous diffusion, namely of the correlated type, charac- 
terized by the following generalized, nonlinear Fokker-Planck equation: 

where (/i, z/) G 7^^, D > is a (dimensionless) diffusion constant, F{x) = —dV{x)/dx is 
a (dimensionless) external force (drift) associated with the potential V{x), and {x,t) is a 
(dimensionless) 1+1 space-time. Let us mention that, in variance with the correlated type we 
are focusing on here, the Levy-like anomalous diffusion is associated with a linear equation, 
though in fractional derivatives 



We intend to consider here a specific (but very common) drift, namely cliaracterized by 
F{x) = ki — {ki E TZ and > 0; A;2 = corresponds to tlie important case of constant 
external force, and ki = corresponds to the so called Uhlenbeck-Ornstein process [p!5[1). 
The particular case fi = u = 1 corresponds to the standard Fokker-Planck equation, i.e., to 
normal diffusion. The particular case F{x) = (no drift) has been considered by Spohn ||^ 
for /i = 1 and arbitrary u (for instance, u = 3 satisfactorily describes a standard solid-on- 
solid model for surface growth), and has been extended by Duxbury [|l^ for arbitrary /i and 
u. The case (/i, fci) = (1,0) has been considered by Plastino and Plastino ^7^. Our present 
discussion recovers all of these as particular instances. 

First, let us illustrate the procedure we intend to follow, by briefly reviewing normal 
diffusion {fj, = u = 1). We wish to optimize 5*1 (given by Eq. (2)) with the constraints 



du p{u) = 1 , (4) 

{u - um)i = J du{u- um) Piu) = (5) 

and 

((m-mm)^)i= I du {u - umY p{u) = a'^ , (6) 



Um and a being fixed finite real quantities. The optimization straightforwardly yields the 
solution 

„-/3(m-ma/)2 

Pi{n) = (7) 

^1 

with 

Z, = Jdu e-^("-"^^)' = (7r//3)i/' (8) 

where (3 = 1/T is the Lagrange parameter associated with the constraint (6) and satisfies 
/3=l/(2a2). 

On the basis of Eq. (7) we propose, for the /z = z/ = 1 particular case of Eq. (3) (i.e., 
the standard Fokker-Planck equation), the ansatz 



with 



m 
m 

It follows straightforwardly that A = 2, 



Zi{t) 



(10) 



m 



and 



hence 



ZijO) 
Zi{t) 



dxMjt) 
dt 



I \ ^1 

^2 



(11) 



ki - k2XM{t) 



(12) 



-k2t 



(13) 



To discuss the ^2 = case we can use e ~ 1 — 2A;2t, which implies xuit) = xm{0) + kit 
and l//3(t) = [l//5(0)] + ADt, which, in the limit t oo (i.e., t » l/[4D/3(0)]), yields the 
familiar result l//?(t) ~ 4Z}t. This result implies, in turn, the celebrated Einstein expression 
{{x — t for Brownian motion. 

Let us now address the general {fi, u) case. Following along the lines of Alemany and 
Zanette |^ (and the generic framework of the generalized thermostatistics [0,0]) we now 
wish to optimize Sg (given by Eq. (1)). The constraints are Eq. (4), 



{U - UM)q = du{u- Um) [P{u)y = 



(14) 



(which generalizes Eq. (5)) and 



{{U - UMf)q 



{u-UMf[p{uW = 0^ 



(15) 



(which generalizes Eq. (6)). This is an appropriate moment for commenting that the reason 
for using [p('u)]'' (instead of the familiar p{u)) in the constraints (14) and (15) is the (very 



essential) fact that by doing so we preserve [|T^] the Legendre structure of Thermodynamics 



and (through the nonnegativity of Cg/q [|1^], where Cg denotes the specific heat) guarantee 
thermodynamic stability. Let us consistently stress that the constraint (14) is equivalent to 
{u)q = {uM)q, but uot to {u)g = um (siuce, unless q = I, {uM)q 7^ um)- All these peculiarities 
are of course originated by the essential nonextensivity that the index q introduces in the 
theory. For example, if we have two independent systems A and B (i.e., Pa*b{ua,ub) = 
'Pa{ua)pb{ub)), we immediately verify that Sq{A*B) = Sq{A) + Sg{B) + {1 — q)Sg{A)Sg{B). 
It is straightforward to see that the above described optimization of Sq yields 

with 

Zg = Jdu[l- f3{l -q){u- UM?]^- (17) 

In the limit g ^ 1, these equations reduce to Eqs. (7) and (8), respectively. The corre- 
sponding ansatz for solving Eq. (3) now is 

with 

A 



m 



(19) 



Zq{t)\ 

(as before, P{t) and Zq{t) are nothing but the scaling of space with time). A tedious (but 
straightforward) calculation yields A = 2/i, and q = 1 + fi — u. An equation for Zg{t) is also 
found, namely 



2uDPiO)[ZqiO)f'^ - k2[Zq{t)r+'' - 

/i d[Zq{t)r+'^ 



0, (20) 



fj, + u dt 

which can be solved by substituting Z(t) = Zg{t)^^^ . The resulting solution is (/or all values 
ofh) 
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with 



K2J K2 



(21) 



K2 = (22) 

2z/D/?(0)[Zq(0)]^-'^ ^ ^ 



and 



T = 



(23) 



A;2(/i + z/) 

(See Fig. 1). The function xuif) is the same as in the case of normal diffusion, Eq. (13), 
since it only describes the motion of the average of the distribution pq{x,t), and does not 
depend on the way in which it spreads. j3{0) and Zq{0) are determined by the initial condition 
(i.e., by pg{x, 0)). 

For k2 = 0, Eq. (21) becomes 

Z,(t) = [lZ,m'--^+'-^^^^^£m?^t\^ (24) 

which, for t 00, asymptotically recovers Duxbury's solution [|l^, namely l/j3{t) oc 
[Zq{t)Y^'- OC t~ . As we see, fi/u = 1, > 1 and < 1 respectively imply that [x{t) — XM{t)]^ 
scales like t {normal diffusion), faster than t (superdiffusion) and slower than t {subdiff'u- 
sion). The limits fi/u = and fi/u = 00 correspond to "no diffusion" and ballistic motion, 
respectively. 

For {fi, ki) = (1, 0), the present set of equations reduces to that of Plastino and Plastino 

Let us mention that the general solution given in Eq. (21) can be derived from that for 
11 = 1 (and arbitrary ki) by defining p{x,t) = [p{x,t)Y, and V = u/fi, as can easily be seen 
from Eq. (3). 

Finally, by using Eq. (19) with A = 2/i, we can verify that 

I dxpq{x,t) = [Zq{t)/Zq{Q)Y-^ I dxpq{x,0) . (25) 
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Consequently, the norm ("total mass") is generically conserved for all times only if // = 1 
(V K2) or if i^^2 = 1 (V/i). For < < 1 (a common case), the norm monotonically 
increases (decreases) with time if > 1 (/x < 1). If > 1, it is the other way around. 

Before ending let us mention that, also when t grows to infinity, the solutions we have 
found must be physically meaningful. This imposes ixjv > —1. Indeed, if k2 ^ 0, t in 
Eq. (21) must be positive, which implies 11/ v > —1. Also, if A;2 = 0, x must scale with an 
increasing function of t] hence, (5{t) must decrease with which implies (through Eqs. (19) 
and (24)) 2///(/x + z/) > 0, hence, the already mentioned restriction applies once again. The 
entire picture which emerges is indicated in Fig. 2 (we have not focused the < region 
because that would force us to discuss the stability of the solutions with respect to small 
departures, and this lies outside of the scope of the present work). 

Summarizing, on general grounds, we have shown that thermostatistics allowing for 
nonextensivity constitute a theoretical framework within which a rather nice unification 
of normal and correlated anomalous diffusions can be achieved. Both types of diffusions 
have been founded, on equal footing, on primary concepts of (appropriately generalized) 
Thermodynamics and Information Theory. 

On specific grounds, we have obtained, for a generic linear force F{x), the physically 
relevant exact (space, time)-dependent solutions of a considerably generalized Fokker-Planck 
equation, namely Eq. (3). 

It is with pleasure that one of us (C.T.) acknowledges warm hospitality by B. Widom 
at the Baker Laboratory. This work was carried out in the research group of B. Widom, 
and was supported by the National Science Foundation and the Cornell University Materials 
Science Center. 
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FIGURES 

FIG. 1. The n/v = 1/3 example: (a) Time dependence of f3{0)/f3{t) = [Zq{t)/Zq{0)]'^t' for 
Zg{0) ^ and typical values of K2 (indicated at the right of each curve). The curve for K2 = lies 
on the vertical axis. For = 0.25, 0.5, and 2 the asymptotic values for f/r ^ 00 are shown by the 
dashed lines, (r is defined in Eq. (23).); (b) Time dependence of {f3{0)[Zq{0)]'^i^}/f3{t) = [Zg{t)]'^i^ 
for Zq{0) = 0, P{0)[Zq{0)]'^^' ^ 0, and typical values of K'2 = k2/{2i^Df3{0)[Zq{0)]'^f'} (indicated 
at the right of each curve). The curve for K2 = 00 coincides with the horizontal axis. All curves 
saturate at a finite value as t — > 00, except that for K'2 = 0, which is proportional to t^i^/ii^+^) for 
all t. 

FIG. 2. "Norm conservation" means that N = J dx pq{x,t) is time-invariant; "Norm creation" 
means that A'" monotonically increases (decreases) with time if < 1 {K2 > 1); "Norm dissipa- 
tion" means that N monotonically decreases (increases) with time if ^2 < 1 {K2 > 1). "Normal 
diffusion", " Superdiffusion" and " Subdiffusion" refer to the fact that, for A;2 = 0, {x — xm)"^ scales 
like t, faster than t and slower than t, respectively. The standard Fokker-Planck equation corre- 
sponds to iJ. = u = q = 1. For the precise meaning of " unphysical" , see the text. On the /x = 1 
line we have q = 2 — v; consequently, when v varies from 00 to —1, q varies from —00 to 3, which 
precisely is the interval within which Eq. (4) (and, consistently, / dx Pq{x, 0) = 1) can be satisfied. 
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